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Abstract 

A part analytical, part numerical ideal MHD analysis of low-frequency Alfven 
wave physics in the H-1 stellarator is given. The three-dimensional, compressible 
ideal spectrum for H-1 is presented and it is found that despite the low /3 of H-1 plas- 
mas (/3 ~ 10~^), significant Alfven-acoustic interactions occur at low frequencies. 
Several quasi-discrete modes are found with the three-dimensional linearised ideal 
MHD eigenmode solver CAS3D, including beta-induced Alfven eigenmode (BAE)- 
type modes in beta-induced gaps. The strongly shaped, low-aspect ratio magnetic 
geometry of H-1 causes CAS3D convergence difficulties requiring the inclusion of 
many Fourier harmonics for the parallel component of the fluid displacement eigen- 
vector even for shear wave motions. The highest beta-induced gap reproduces large 
parts of the observed configurational frequency dependencies in the presence of hollow 
temperature profiles. 

1 Introduction 

A rich low-frequency wave phenomenology showing Alfvenic characteristics has been ob- 
served in plasmas generated on the H-1 heliac.[Tl [21 E] Previous ideal MHD studies of 
these waves in cylindrical geometry has shown that cylindrical shear Alfven continuum 
turning points reproduce important features of measured magnetic frequencies [U |2], and 
cylindrical global Alfven eigenmode (GAE) modelling produces plausible radial eigenmode 
profiles and magnetic attenuation predictions [3]. However, it is well known that Alfven 
waves existing independently in a cylindrical plasma will couple in a toroid due to mag- 
netic field strength modulations, and the same is true for coupling between acoustic waves 
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and Alfven waves at low- frequencies. [1] Since these coupling effects are expected to feature 
prominently in low aspect-ratio, strongly-shaped plasmas like those created on H-1, this 
paper aims to provide a geometrically accurate treatment of H-1 wave physics within the 
theoretical framework of ideal MHD. 

Two significant drivers of Alfven wave activity studies in the magnetic confinement 
fusion context are the direct effects such wave activity may have on plasma confinement, 
through energetic particle interactions for example, and the diagnostic utility wave activ- 
ity provides as proxy for determining underlying plasma parameters. [1] In particular, low- 
frequency Alfven-acoustic waves have drawn increased attention in recent years as they 
are generally poorly understood. [HI E] Even under the simplifying assumptions of ideal 
MHD, the low-frequency domain is difficult to treat theoretically because, in the presence 
of geodesic curvature, couphng between Alfven and sound branches introduces considerable 
complexity to the linearised eigenfrequency spectrum. Furthermore, ideal MHD provides 
a poor description of compressible motions|7], and it omits thermal ion interactions, such 
as ion Landau damping, which become significant in this range of frequencies. [8J How- 
ever, ideal MHD has the important advantage that it is simple enough to allow for an 
accurate treatment of the geometry, providing a good first model for understanding the 
low- frequency domain. [H [9l [10] 

In this work we present the three-dimensional, compressible ideal spectrum for H- 
1, based on numerical simulations with the three-dimensional continuous spectrum code 
CONTI[lT] in combination with the three-dimensional ideal MHD linear eigenmode solver 
CAS3D.[12l [13] Significant coupling of eigenmode Fourier harmonics in both toroidal and 
poloidal directions are induced by a substantial mirror term (inherent in a helical axis, low 
aspect-ratio configuration) in combination with a variety of strong field-strength modula- 
tions within a cross-section. We have found that the conditions for convergence of CAS3D 
in Fourier space are quite stringent for H-1 configurations, in the sense that shear Alfven 
frequencies dip sharply near the core unless the component of fiuid motion along the field 
lines has a Fourier basis including at least two toroidal side-bands, even for high-frequency 
modes with negligible acoustic interaction. An explanation for this finding is given. 

The first two helical Alfven eigenmode (HAE) gaps and the toroidal Alfven eigenmode 
(TAE) gap are identified. Significant beta-induced gap structure is found, starting at 
around 35 kHz. Consequently, H-1 fiuctuations, which typically lie near or below 35 kHz, 
may be considered sub-GAM modes. [5] Based on analytic estimates, it is proposed that 
H-1 fiuctuations may include beta-induced Alfvenic eigenmodes partially reproducing char- 
acteristic "whale-tail" structures in configuration space due to hollow temperature profiles. 
Low frequency discrete modes are presented. 

The outline of the paper is as follows. Section 2 provides an overview low-frequency 
H-1 wave activity. Section 3 discusses H-1 equilibrium and magnetic geometry. Section 4 
gives a summary of low-frequency stellarator Alfven wave physics as well as a discussion of 
CAS3D convergence requirements for H-1. Section 5 presents the ideal spectrum and low- 
frequency discrete modes. In Section 6 we show that, in the presence of hollow temperature 
profiles, the highest beta-induced Alfven eigenmode (BAE) gap frequency reproduces parts 
of the configuration-space whale-tails. We conclude in Section 7. 
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2 Overview of low- frequency waves in H-1 plasmas 



In this section we discuss low-frequency wave activity on H-1 and give a brief overview 
of previous studies. H-1 is a "flexible" lieliac, which means that a variety of magnetic 
geometries can be accessed by adjusting Kh, a parameter measuring relative coil currents 
(held constant during a discharge), where < k/^ < 1.2. Different values of correspond 
to different plasma parameters, flux surface geometries and rotational transform h = 1/q 
profiles (see ref. ^ for details). Perhaps the most important of these is the * dependence 
on Kh, shown in Figure [T} 
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Figure 1: : Rotational transform a as a function of average minor radius for selected 
values of k^, where the plasma edge lies at r = 0.2 m. Low-order rational values are 
indicated on the right. 

The configurational dependence of fluctuation frequencies is our main focus as H-1 
discharges are usually dominated by a single coherent steady-frequency fluctuation present 
in a time interval with limited variation in equilibrium magnetic field strength and plasma 
density. PQ |3] Figure |2] shows the frequencies measured by poloidal Mirnov coils at different 
times within discharges for a selection of discharges. [2] 

It has been shownlU [2] that the "V" structures (which we term "whale-tails") near 
Kh = 0.4 and Kh = 0.75 in figure [2] qualitatively follow the turning points of the cylindrical 
shear Alfven frequency oja = \n — *m|i?o/-R/ioP for ijn^n) = (4,5) and {m,n) = (3,4) 
respectively (compare with low-order rationals in Figure [T]), where the turning points of 
ua depend on Kh through a and p (the mass density), with the average magnetic field 
strength Bq ^ 0.46 T held constant during a discharge. Mode-number assignments for 
these two regions are obtained from poloidal Mirnov array phase signatures [H [2] and initial 
analysis of helical Mirnov array measurements, where m and n are the azimuthal and axial 




3 



100 



80 - 




J" 



I 




0.2 



0.4 



0.6 



0.8 



1.0 



1.2 



Figure 2: : Magnetic field fluctuation frequencies measured by poloidal Mirnov array 
coils against Kh. 

mode-numbers respectively. These are dominant mode-number assignments as multiple 
mode-numbers have not yet been resolved from the available data. Assuming the best 
information available for mass density and using H-l's average major radius R = 1 m 
in the formula for ua, predicted frequencies were found to be a factor of A ~ 3.7 times 
larger than measured frequencies. The discrepancy factor is reduced to A ~ 3 if the helical 
extension of the arclength of the magnetic axis is accounted for by setting R = 1.25 m. 
This discrepancy was not the same for each of the four branches of the two whale-tails 
and mode localisation assumptions were made for some configurations. Consequently, the 
value A ~ 3 should be taken as a coarse measure of the theory-measurement frequency 
discrepancy. 

The frequency scaling behaviour, which is qualitatively analogous to frequency sweep- 
ing of reversed-shear Alfven eigenmodes (RSAE) as q changes in tokamak discharges [Ij, 
suggests that the low-frequency modes are CAEs and non-conventional CAEs (NGAEs) at 
local minima and maxima in ua respectively. The low toroidal currents of H-1 favour the 
existence of NGAEs, whereas RSAEs have only been observed in tokamaks during "Alfven 
cascades" in which several waves with different mode numbers sweep through a broad fre- 
quency range during the course of a discharge, unlike the steady, solitary H-1 modes. [H] 
However, no robust explanations have been found for the quantitative discrepancy. While 
significant impurities are likely to be present in H-1 plasmas, cylindrical frequencies are 
expected to depend weakly on this effect as oja is proportional to the inverse square root of 
the effective mass. Collisions with neutrals are likewise not thought to be able to explain 
such a large discrepancy. pLj 

Cylindrical ideal MHD modelling of CAE eigenmodes under H-1 conditions has shown 
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agreement with magnetic attenuation and optical emissions measurements, p] However, 
as that study assumed cyhndrical geometry, no comparisons were made with modes other 
than GAEs. Furthermore, that study highhghted the multiply-peaked nature of the Fourier 
spectrum of magnetic phase measurements, pointing to the importance of mode-number 
coupling in H-1. We will return to this "GAE hypothesis" in view of the 3-dimensional 
spectrum. 

3 H-1 equilibrium with VMEC 

We use a Boozer coordinate representation (s, () for ideal MHD equilibrium. Here 
s denotes the normalised toroidal flux ip/ipedge where 27r'?/' is the toroidal flux, and d,C, 
are magnetic poloidal and toroidal Boozer angles respectively. The Variational Moments 
Equilibrium Code[Tn| (VMEC) solves the ideal equilibrium equation = (V x B) x 
B in 3-dimensions assuming that nested flux surfaces fill the plasma volume. VMEC 
chooses the poloidal angle 6'vmec so that the spectral width of the Fourier representation 
6'vMEC5 <P)i Z{s, ^vMEC, for each magnetic surface is minimised [16] , where {R, 0, Z) 
are cylindrical coordinates with the geometric toroidal angle. This yields significant com- 
putational benefits, although further computation is required to obtain the larger magnetic 
coordinate Fourier representation {R{s, {}, (), Z{s, {}, ()). 

We have opted to give VMEC predetermined t profiles since the profiles shown in 
Figure [1} which are vacuum field profiles, are considered reliable when a plasma is present 
due to the fact that H-1 operates at low-/3 (/3 ~ 10~^).[17j Pressure profiles are obtained 
from electron density measurements using the ideal gas law p = UekBTe + riikBTi = 2nkBTe 
assuming = Ue and constant plasma temperatures Ti = = 20 eV. |2] We have run 
VMEC in fixed-boundary mode, with boundary surfaces obtained using DESCUR which 
generates the representation {R{s,9yuec,(P), Z{s,9ymec,(P)) from a field-line trace on a 
magnetic surface. [16j The field-line trace has field periodicity of = 3 (i.e. all equilibrium 
quantities are 3-periodic in the toroidal direction) and is stellarator symmetric 

R{s, = R{s, -0 Z{s, = -Z{s, -0 (1) 

as symmetry-breaking terms arising from the response of the magnetic field coils to the 
applied field are ignored. 

Transformation to Boozer coordinates is performed with the MC3D code (part of the 
CAS3D package). This transformation gives poorly reconstructed pressure profiles (from 
Vp = (V X B) X B where B is reconstructed in the Boozer coordinate Fourier basis) for 
H-1 configurations unless a VMEC equilibrium with a large number magnetic surfaces Ng 
is supplied. Figure |3^ shows the improvement of the agreement between the prescribed 
VMEC pressure gradient and the MC3D-reconstructed pressure gradient with increasing 
Ng. The reconstructed pressure for the Ng = 500 case is shown in figure |3|3 along with 
the prescribed VMEC pressure and electron density inversions from a Kh = 0.54 discharge. 
Singular behaviour at the core is a consequence of the s = coordinate singularity arising 
from the vanishing of the poloidal flux gradient, which affects the VMEC equilibrium as 
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well as the Boozer coordinate transformation. [15] The Vp discrepancy near the edge may 
be due to the difficulty of representing B on the more strongly shaped outer magnetic 
surfaces, even though the geometry of the flux surfaces is accurately reconstructed. The 
largest A^^ = 500 Fourier table used before memory requirements became prohibitive was 
< m < 40 and —30 < n < 30. Figure |3^ shows little improvement moving from 
Ns = 350 to Ns = 500, although agreement is harder to obtain for higher /3 values (for 
example /3 = 5 x 10~^ compared with the present /3 = 1 x 10~^) where Ng = 500 flux 
surfaces or more are necessary (larger /3 values are expected with the recently upgraded 
H-1 heating system [IB]). 



(a) (b) 
Vp P 




Figure 3: : (a) Prescribed VMEC pressure gradient compared with the MC3D- 
reconstructed pressure gradient for increasing numbers of flux surfaces A^^ = 
100,200,350,500. (b) Pressure reconstructed from the A^^ = 500 case in (a), compared 
with the prescribed VMEC pressure as fitted to electron density inversions [19j. All pres- 
sures normalised to core fitted pressure. This is a Hh = 0.54 configuration. 



Stellarator symmetry simplifies the Fourier representation of equilibrium quantities by 
permitting sine-only or cosine-only expansions. The Boozer coordinate Fourier expansion 
oi B = |B| on a field period can be written as [151 IE] 

B{s, ^X) = Y1 ^™"(^) cos(m79 + nNQ = 5o 1 1 + 6^„(s) cos(mt9 + nNQ j • (2) 

mn \ mn / 

Figure |4^ gives the flux surface dependence of the B^n for = 0.54 in magnetic co- 
ordinates showing that m = 1 Fourier components dominate the poloidal field strength 
modulation. Similarly, the Boozer coordinate Fourier representation of the contravariant 
metric component 5''^* = Vs ■ Vs can be written as 



g^%s, ^,0 = Y1 a'r^nis) cos(m^ + nN() = | ^ + E 9mn{s) cos(m^ + riNQ ] . 

mn 'Vedge 



(3) 
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where the last expression on the right defines gmni^) = 5o('0) which is related to 

the elongation of flux surface cross-sections. [E] The flux surface dependence of the g^^^ are 
shown in Figure |4}d, illustrating the strong ellipticity (m = 2) and triangularity (m = 3) 
of the bean-shaped cross-sections. 




Figure 4: : (a) Fourier components of the magnetic fleld strength B = |B| in Boozer 
coordinates on a fleld period. The largest four components shown with thick lines are: 
(m,n) = (0, 0)-red (with Bq = 0.46 subtracted), (l,0)-blue, (1, — l)-green, (1, — 2)-orange. 
All other components in thin black, (b) Fourier components of the contravariant metric 
component g'^^ = Vs ■ Vs in magnetic coordinates. The largest three components shown 
with thick lines are: {m,n) = (0,0)-red, (2, — 2)-yellow, (3, — 3)-brown. All other compo- 
nents in thin black. This is a. Kh = 0.54 conflguration. 



4 Wave Theory Preliminaries 
4.1 Mode Families 

The toroidal A^-periodicity of stellarators leads to the notion of independent mode-families. 
A choice of mode family is the stellarator equivalent to a choice of toroidal mode number 
n for tokamaks, as different n do not couple in axisymmetric geometry. 

The condition for coupling between Fourier mode-numbers can be obtained from a vari- 
ational formulation of the linearised ideal MHD equations. [20] Eigenvalues u and plasma 
displacement eigenfunctions ^ of the linearised ideal MHD equations extremize the La- 
grangian 

C = u'^K - 5W (4) 
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Table 1: 
□ . 



Summary of H-1 mode-families with = represented by ■ and iV = 1 by 



with respect to variations of ^|20], where the kinetic energy is K 
potential energy 5W can be written as 



5W 



\ I Pl^l'^dr and the 
C|2-Z}(^ ■ Vs)2 + 7p(V-0'f^r, (5) 



where 



D 



Vx(^xB) + :^^^V., 

I Vs|"' 

,(J X Vs) ■ (B- V)Vs 



(6) 
(7) 



and boundary and vacuum energy terms have been neglected i.e. the plasma is assumed to 
be surrounded by a perfectly conducting wall. The latter assumption is made for simphcity, 
as only fixed boundary equilibria are considered here, and it does not affect the mode family 
argument that follows. 

All terms in the 6W integrand are of the form "perturbation x perturbation x equilib- 
rium" . Substituting the general (combined cosine and sine) Fourier expansion for ^ (with 
poloidal mode-numbers m and toroidal mode-numbers n) and the general Fourier expan- 
sion for each equilibrium quantity (where the toroidal mode-numbers have the form nN), 
the orthogonality of the trigonometric functions leads to the condition for coupling 



n + n' E N7j or 



n 



n' E NZ, 



where NZ = {. . . , -2A^, -A^, 0, iV, 2A^, . . .}.[I2] 

The N = 3 toroidal periodicity of H-1 configurations implies the existence of two 
independent mode-families = 0, 1, shown in Table [T| which are uncoupled. We will 
focus on the N = 1 mode family since it contains both of the dominant mode-number 
pairs (4,5) and (3,4) for the two whale-tails (see section [2]). 

4.2 Continuous Spectrum 

The continuous ideal MHD spectrum for a general three-dimensional toroid can be obtained 
by CONTI[II] which solves the coupled continuum equations [2]j| 



^,(B.V,(^(B.V) 



IP + B' IP 1 



e. + 7P«:s(V-^) =0, 

(v-0 = o, 



(9) 
(10) 



8 



where 7 is the adiabatic index, uj is the angular frequency, = ^ ■ (B x V'?/')/! VV'p is 
the surface displacement and (B X ViIj)/B'^ is the geodesic curvature, which is 

proportional to the surface component {Vip xh) ■ k/B of the magnetic field line curvature 
K = (h ■ Vb) where b = B/B. The continuum equations are derived from the linearised 
ideal MHD equations. For a given flux surface ip they give the condition for the existence 
of non-square-integrable eigenfunctions with a spatial singularity at ip (see ref. [21] for 
details and further references). 

If Kg vanishes, equations ^ and (10) decouple and describe the shear Alfven and 



compressional (acoustic and compressional Alfven) branches of the ideal MHD continuum 
respectively. Coupling between Fourier components of the shear Alfven spectrum arises 
due to modulation in iVipl/B, where Fig. |4] shows the numerator (up to a constant factor) 
and denominator in this expression. 

In the presence of nonzero p and k^, equations ^ and (10) couple. Setting aside the 



high-frequency compressional Alfven waves, it can be seen from the identity 

V-^ = b-Veii+V-^^ (11) 

that this coupling involves two factors. The first term on the right represents coupling 
between the acoustic motion parallel to the magnetic field lines and the surface displace- 
ment ^s- The second term represents the compressional response due to motion across the 
equilibrium field lines = B x E/i?^ in the presence of geodesic curvature. jH] This effect 
is closely related to the electrostatic geodesic acoustic mode (GAM). [23] 

Due to the mixing of shear and compressional motion at low frequencies the clean 
division between them breaks down and it is necessary to introduce classification schemes. 
Three schemes are used in this paper, classifying eigenfunctions as shear Alfvenic if: (1) 
the largest Fourier harmonic of V ■ ^ is larger than the largest Fourier harmonic for C,s in 
Boozer coordinates (CONTI classification), (2) the parallel component K\i = ^ J p\^\\\'^dr of 
the displacement kinetic energy K dominates i.e. K\\/K > 0.5 (polarisation classification) 
and (3) the plasma compression energy ^ f 7p(V ■ $,)'^dr dominates the other terms (field 
line bending and compression energies) in SW (potential energy classification). 

A criterion for estimating the strength of the Alfven-acoustic interaction has been 
proposed, given in terms of a "sound parameter" , [H] 

® = y-^, (12) 

where e = ^2ijj/Bo/R is a measure of the inverse aspect ratio (since the toroidal flux is 
27rilj = BoTrr'^ in a torus) and 610 and 60 > 1 are defined in equations ^ and ^ respectively. 
If C5 ^ 1, the interaction is predicted to be weak and the Alfven frequency is zero at the 
surface where n— hm = as expected for the incompressible shear Alfven spectrum. On the 
other hand, if -C 1, significant deviations from the incompressible shear Alfven spectrum 
are expected with the formation of Alfven-sound gaps and in some cases, immediately above 
these, /5-induced gaps bounded above by geodesic acoustic frequencies ug located near the 
surface where n — *m = 0.[l3] For H-l the sound parameter ranges from 35 ~ 0.75 at 
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Kh = 0.3 to ~ 1.1 at Kfi = 0.9 suggesting non-trivial Alfven-acoustic interactions at low 
frequencies. H-1 is unusual in this regard, since e/bio > 1 for most stellarators (H-1 has 
e/^io < 1)) and therefore it is often the case that 3> 1 for stellarators with * > l.[T^ 

An approximate expression for the Alfven resonance ignoring Alfven-Alfven interactions 
(equation (48) in ref. [H]) can be used to obtain estimates for ug by taking the limit 
n — im — > yielding for H-1 

,2 _ x2^2/p2 + ,.,2 _ x^2^2/p2 + ,.,2 _ ^0/\r_ x^2^2/p2 R„^2 ' ^^'^) 



00^ - i^cl/R^ uj^-{N- tfcl/R'' uj^ - {2N - i)^cl/R^ Bod, 



where all flux-surface quantities are evaluated at the surface where n — tm = 0, N = 3, 
Cs = \/lvl P is the adiabatic sound speed, (5o ~ 1-5 and only the largest harmonics identified 
in Fig. |4]with thick lines are considered. The geodesic frequencies ujc form a subset of the 



solutions bj to equation (13). [H] 



4.3 Discrete Spectrum and CAS3D convergence 

The full ideal MHD linearised eigenmode spectrum, including discrete or quasi-discrete 
modes existing within gaps in the continuous spectrum, can be found with CAS3D [T2| [13] 



which solves the eigenvalue problem (5£ = (see section 4.1). 

For H-1 configurations, CAS3D produces poorly converged spectra with shear Alfven 
frequencies dipping sharply near the core unless the parallel displacement is represented 
by many Fourier harmonics. Some insight into this problem can be obtained by rearranging 



equation (11) and solving in magnetic coordinates to give 



n — hm 

where the Jacobian is = [V('?/') • {Wd x VC)]~^- For a physical divergence V ■ ^ that 
needs to be numerically approximated, it can be seen that the spectral width of must 
be larger than that of the other displacement components due to the modulation in B^. 
CAS3D uses Boozer coordinates for which B^ oc 1/5. 

Inadequate representation of has deleterious effects well outside of the frequency 
domain in which parallel ion-acoustic motion is physically important. For simplicity, let us 
consider the case of approximately "pure" shear Alfven waves which should have V ■ ^ ~ 
and ~ 0. The parallel displacement appears only in the fluid compression term 
7p|V ■ and the parallel component of the kinetic energy i^y = | J pl^ypdr. If is not 
represented with a sufficient number of Fourier harmonics, both V ■ ^ ~ and ^y ~ will 
not be well satisfied, and both j 7p(V ■ ^)^(ir and K\\ will be larger than they should be. 
Although 7p|V-^p will still be small relative to the magnetic terms in 5W when /3 is small, 
an infiated K\\ will cause significant distortions to the shear Alfven spectrum. 

This effect has been documented for unstable modes with the 3-dimensional ideal 
spectral code TERPSICHORE. [21] For CAS3D with W7-X equilibria it was found that 
a one-third larger Fourier table for the parallel displacement was sufficient to avoid these 
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problems. [T3] For H-1, due to a significant mirror-term modulation of around 20%, CAS3D 
Alfven spectra are grossly distorted near the core where poloidal couplings are small and 
V ■ ^ ~ cannot locally be satisfied with the inadequate "toroidal representation" of ^y, 
unless the first and second toroidal side-bands n — 2N, n — N, n + N and n + 2N of 
near-resonance harmonics n are included for ^y. The inclusions of these side-bands and 
their poloidal neighbours leads to Fourier tables twice or three times as large as those for 
the other components. Our simulations with H-l-like reduced-mirror-term configurations 
which retain strong poloidal modulations has shown that only the first toroidal side-band 
is required for ^y. 

5 Numerical Results 

This section starts with a description of the broad features of the H-1 linearised spectrum 
before moving on to a more detailed study of the low-frequency part of the spectrum. 
Configuration space can be divided into four regions of interest corresponding to the left 
and right sides of the two whale-tails. On the left side of the Iow-k^ whale-tail, the 
(m, n) = (4, 5) shear Alfven continuum intersects the low-frequency domain because there 
is a flux surface in the plasma with a = 5/4 whereas on the right side of the \ow-Kh whale- 
tail there is no surface for which a = 5/4 and the {771,11) = (4,5) shear Alfven continuum 
rapidly moves to high frequencies with increasing Kh- The same pattern repeats itself for 
the (m, n) = (3, 4) shear Alfven continuum and the high-K/i whale-tail. Consequently, 
we focus our attention on = 0.30 as representative of low-frequency Alfven-acoustic 
interactions on H-1. 

5.1 Alfven spectrum 

The shear Alfven gap structure for a. Kh = 0.30 H-1 configuration obtained with CONTI 
can be seen in Figure [5^. Two HAE gaps are present, a large one at around 150 kHz 
corresponding to a coupling displacement of (Am, An) = (2, 3) (i.e. coupling between 
(m,n) and (m -|- Am, 71 + An)), and above it, a small gap corresponding to (Am, An) = 
(3,3). Above both of these gaps lies the TAE gap at around 350 kHz. This overall gap 
structure is found to be qualitatively the same for all kh, with significant quantitative 
differences when there are large changes in p near Kh = 0.4 and Hh = 0.75. [T] These gaps 
are far too high to offer an explanation for the observed low-frequency wave activity. 

Figure [5|d shows a comparison of the (4, 5) and (3, 4) Alfven branches with their cylindri- 
cal equivalents for Hh = 0.54 (at the overlap of the two whale-tails). This shows that apart 
from the appearance of gap structure, at frequencies greater than 50 kHz the 3-dimensional 
shear Alfven frequency turning points are about a factor of 0.85 lower than correspond- 
ing cylindrical frequency turning points, implying a reduction of the CAE-measurement 
frequency discrepancy to A ~ 2.5 (see section [2]). Ceometric effects therefore alter CAE 
frequencies more than previously thought [H |2], but are not sufficient to bring them in line 
with measurements. 
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(a) 




Figure 5: : (a) CONTI continuous spectrum for N = 1, = 0.30 with with p and p 
hnear in s and constant temperature. Sound branches in grey, Alfven branches in colour, 
(b) Continuous spectrum for = 1, = 0.54 with p and p quadratic in s as in Fig. 
[3jD and constant temperature. Modes classified as acoustic not shown, Alfven branches 
in yellow apart from (4, 5)-red and (3, 4)-blue with dashed lines showing corresponding 
cylindrical frequencies. For both figures physical parameters are B = 0.46, core mass 
density po = 1 x 10~^ kg/m^, core pressure po = 16 Pa and CONTI numerical parameters 
are 400 flux surfaces, 1 < m < 18 and —2 < n < 23 (with multiples of 3 removed to 
exclude N=0). 



5.2 Alfven-acoustic spectrum 

We now turn our attention to the low- frequency domain. Figure |6^ shows the low-frequency 
detail of figure M where the (m, n) = (4, 5) branch has been identified, confirming the 



significant Alfven-acoustic interactions predicted in subsection |4.2[ Several Alfven-acoustic 
gaps are formed, with large ^-induced gaps appearing above some of these notably at 9 kHz, 
12 kHz. Another smaller /3-induced gap is found at 30 kHz which is more easily seen in 



Figure [6p. The positive solutions to equation (13) at s = 0.55 are u ~ 9.4, 13.2,31.1 kHz 
which are therefore all seen to be approximate geodesic frequency solutions. 

Figure [7^ shows the full Kh = 0.3 low-frequency ideal MHD spectrum produced by 
CAS3D, where a small mode-number table (adjacent poloidal mode-numbers and two 
toroidal sidebands) has been used to make the gap structure clear. The CAS3D spec- 
trum differs from the CONTI spectrum in that it shows a strong geodesic cut-off near 
35 kHz and the lower sections of the shear Alfven spectrum are absent. 
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(a) 




Figure 6: : (a) Low frequency detail of figure [5| where the (m, n) = (4, 5) shear Alfven 
branch is coloured red. (b) Detail of (a). 

Several "quasi-discrete" modes are found below 35 kHz, meaning that the mode spatial 
structures have large sections apparently free from continuum interactions. Figure [8^ shows 
an example BAE using a small Fourier basis in order to suppress acoustic resonances. 
Even with the inclusion of many sound branches, eigenmode clusters with shear Alfven 
characteristics are found as shown in [7jD using the potential energy classification scheme. 
Quasi-discrete modes are found in some of these clusters, such as the mode shown in figure 
[SjD, found near the cluster circled in[7)D. As these are fixed-boundary simulations, the fluid 
displacement orthogonal to the flux surfaces must vanish at the edge, which it does due 
to a sharp gradient discontinuity near the edge presumably due to continuum interactions. 
Modes with a very similar structure but the discontinuity appearing further inside the 
plasma are also found. The modes shown in figure [8] may therefore experience significant 
continuum damping. 

Low-frequency quasi-discrete modes dominated by (m, n) = (3, 4) are also found for 
the left branch of the high-Kh whale-tail. Two n = 0.60 quasi-discrete modes are shown in 
figures |9^ and|9]D, with frequencies of 38 kHz and 34 kHz respectively. The mode in figure 
[9]d shows clearly a typical example of the aforementioned near-edge resonances. On the 
right sides of the whale-tails, near k = 0.50 and k = 0.90, we have not been able to find 
low-frequency Alfvenic discrete modes dominated by (m,n) = (4,5) and (m,n) = (3,4) 
respectively. 
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(a) 




Figure 7: (a) Low frequency CAS3D spectrum where the (m, n) = (4, 5) shear Alfven 
branch is coloured red (polarisation classification). Only adjacent poloidal harmonics and 
two toroidal sidebands of the (4, 5) harmonic included in the Fourier table for all mode 
vector components (total of 15 harmonics for each), (b) Same as (a) with larger Fourier 
table and potential energy classification. The components have 16 near-resonance 
Fourier harmonics with m^[^ = 3, mmax = 5, and n^[^ = —10 and ^max = while the 
parallel displacement has mmin = 2, mmax = 6, nmin = — 14 and n^ax = 4. In both figures 
1000 fiux surfaces are used, and Kh = 0.30. 



6 Alfven-acoustic whale-tails? 

Since /3-induced gaps below 50 kHz and quasi-discrete modes with reduced continuum 
damping exist in these gaps, we offer a partial alternative explanation of H-1 wave activity 
in terms of /3-induced Alfven eigenmodes (BAEs) or /3-induced Alfven-acoustic eigenmodes 
(BAAEs). 



Figure 10 a, shows the dependence of the (m, n) = (4, 5) cylindrical shear Alfven 



continuum ua near k = 0.4 (the low-/? whale-tail) and figure 10 d shows the same for 
{m,n) = (3,4) near k = 0.75 (the high-K whale-tail). The left branches of both of 
these whale-tails coincide with the inward migration of the outermost surface at which 
the cylindrical shear Alfven continuum vanishes Sres- The geodesic acoustic minimum of 
the 3-dimensional compressible shear Alfven spectrum is located near Sres (compare figure 
|6|i with the k = 0.30 curve in figure ^Ok), and it will migrate inwards with Sres- 

There are indications that H-1 electron temperature profiles increase steeply towards the 
plasma edge and also that they "hollow out" as = 0.4 and Kg = 0.75 are approached. [19] 
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Figure 8: (a) CAS3D discrete mode circled in Figure [7^. (b) CAS3D discrete mode 
found near cluster circled in Figure [TJd. Both cases are for k = 0.30. Some mode harmonics 
have been labelled and it is seen that both modes are dominated by {m,n) = (4,5). The 
vertical axis is the normalised fluid displacement orthogonal to the flux surfaces 



Consider the left side of the \ow-Kh whale-tail covering the range 0.3 < k/i < 0.4. In the 
absence of ion temperature measurements or accurate spatially-resolved electron tempera- 
ture profiles for these values of Kh, we postulate a total temperature profile T(k/i, Srcs) that 
evolves linearly from T = 27 eV at = 0.3 and s^^s ~ 0.55 to T = 8 eV at = 0.38 
and Sres ~ 0.25 where the values Sres are taken from the {m,n) = (4,5) Alfven branch 
(compare with figure 10). We emphasise; this temperature profile is not intended to rep- 
resent the temperature for a single discharge, it represents local temperatures in (s, Kh) 
space postulating a simple linear profile consistent with preliminary spatially resolved elec- 
tron temperature reconstructions [12], assuming = These temperatures can be used 
in equation (13) to estimate how the ~ 32 kHz branch of evolves with k/j for the 
{m,n) = (4,5) shear Alfven branch in the range 0.3 < /t/i < 0.4, which is shown in figure 
[TTj Comparison of figure VL with figure |2] shows that BAEs or BAAEs excited just below 
ug will show a very similar configurational frequency dependence with the left side of the 
low-K whale-tail. 

Unfortunately the available electron temperature measurements do not extend to the 
outer flux surfaces (s > 0.6) which would be necessary for an explicit comparison of ug in 
the range 0.6 < Kh < 0.75 with the left side of the high-K whale-tail. It is clear though 
that a similar frequency dependence on will be observed provided that the temperature 
increases accordingly; that it increases to some degree is fairly certain. 

The quasi-discrete modes found in section 5^ will probably experience enhanced con- 
tinuum interactions in the presence of the non-constant temperatures used in this section. 
Whether these interactions will outweigh the relevant driving mechanisms is an important 
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Figure 9: The normalised fluid displacement orthogonal to the flux surfaces for: (a) 
A {m,n) = (3,4) discrete mode at 38 kHz and (b) a (m, n) = (3,4) quasi-discrete mode at 
34 kHz. Both cases are for k, = 0.60. 

question that requires more detailed physics than ideal-MHD alone provides. [lO] 

7 Discussion and Conclusions 

The linearised ideal MHD spectrum is the starting point for analysis of suspected Alfven 
wave observations, as the gap structure of that spectrum identifles frequency ranges with 
reduced continuum damping. We have presented the three-dimensional, compressible ideal 
spectrum for H-1 plasmas. The presence of /3-induced gaps and partial gaps, the existence 
of discrete modes in these gaps and the predicted conflgurational dependence of the highest 
BAE gap frequency in the presence of H-l's hollow temperature proflles can account for 
important features of wave frequency measurements, namely the left sides of the V-shaped 
"whale-tail" frequency structures observed in conflguration space. Thus, a potentially 
complete explanation of the left sides of both whale-tails has been provided. 

This "BAE hypothesis" does not, however, provide an explanation for the right sides 
of the whale-tails as the relevant resonant surface will leave the plasma and the Alfven 
continuum will rapidly approach higher frequencies. For example, at Kh = 0.54 the (4, 5) 
minimum frequency (see Fig. |5}d) is more than twice as large as measurements. Further- 
more, no discrete modes were found below the (4,5) and (3,4) shear Alfven continuum 
minima near k = 0.45 and k = 0.9 (the right-hand sides of the whale-tails). Nevertheless, 
as quasi-discrete Alfvenic modes with dominant mode numbers (4, 5) and (3, 4) could not 
be found at low-frequencies for k = 0.45 and n = 0.9 respectively, there do not appear to 
be any ideal MHD candidates for the right sides of the whale-tails apart from GAEs that 
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Figure 10: (a) Cylindrical shear Alfven continua for (m, n) = (4, 5) showing Hh dependence 
near = 0.40. (b) Cylindrical shear Alfven continua for {m,n) = (3,4) showing 
dependence near Kh = 0.75. 

have had their frequencies lowered by non-ideal effects, or by an effective mass rrics around 
six times greater near the plasma core than our impurity-free assumption of meg = 2.5mp 
where rrip is the proton mass. The factor of six follows from the GAE-measurement dis- 
crepancy factor which we showed was reduced to A ~ 2.5 when 3-dimensional effects are 
included. 

CAE and BAE wave physics should be experimentally distinguishable by the fact that 
only BAEs are strongly affected by plasma temperature. This highlights the importance 
of rehable spatially-resolved temperature measurements which are the subject of future 
work. It is important that uncertainty in the effective mass and its radial dependence be 
reduced, as this parameter affects both shear Alfven and acoustic frequency predictions. It 
is hoped that measurements from the new helical Mirnov array have the potential to provide 
multiple mode-number assignments as well as eventually allowing for comparison with 
eigenmode polarisation predictions. It may also be possible to obtain accurate spatially 
resolved measurements of the relative plasma density fluctuations which, when compared 
with Mirnov signals would help in characterising shear Alfven versus acoustic components 
of the observed wave dynamics. Ultimately, candidate modes should be compared with 
Mirnov array signals and internal measurements for a conclusive identification. 

Efforts are also underway to obtain energetic particle drive estimates for CAS3D eigen- 
functions and extend equilibrium and spectral calculations to incorporate a free boundary 
and vacuum field. Ideal MHD must be interpreted with due caution at low frequencies due 
to the poor description it provides of compressible motions[7j, and the omission of thermal 
ion interactions, such as ion Landau damping. For these reasons it may prove fruitful to 
pursue kinetic models of sub-GAM H-1 wave physics. O [8] 
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Figure 11: The ~ 32 kHz branch of ojg as a function of Hh based on a postulated hollow 
total temperature profile. 
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